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Using chiral anomaly, we discuss the pair creation of massless fermions in an electric flux tube E 
under homogeneous magnetic field B parallel to E. The tube is axial symmetric and infinitely long 
in longitudinal direction. In the limit B ^> E, we can analytically obtain the spatial and temporal 
behaviors of the electric field and azimuthal magnetic field generated by the produced fermions. We 
find that the life time t c of the electric field is shorter as the width of the tube is narrower. Applying 
it to the glasma in high-energy heavy-ion collisions, we find that color electric field decays fast such 
as t c ~ QJ 1 with saturation momentum Q s . 

PACS numbers: 12.38.-t, 24.85. +p, 12.38.Mh, 25.75.-q 
Schwinger mechanism, Chiral Anomaly, Color Glass Condensate 

Pair creation of charged particles in a classical strong electric field has been theoretically discussed[l[ for long time 
since the discovery of Klain paradox, although such phenomena have not yet been observed. The particle creation is 
known as Schwinger mechanism 2-]. Here, we are addressed with the problem in a quite different way from previous 
ones and apply our results to color gauge fields (glasma) produced in high-energy heavy-ion collisions. 

The glasma have recently received much attention. The gauge fields are classical longitudinal color electric and 
magnetic fields. It is expected that the decay of the glasma leads to thermalized quark gluon plasma (QGP) observed 
in recent RHIC experiments. It is a possible mechanism that the color electric field decays with the pair creation of 
quarks. Indeed, the color electric field has been phenomenologically assumed in previous papers [3| to be generated 
in the collisions and to lead to QGP owing to the pair creation of quarks. On the other hand, the generation of the 
classical gauge fields has recently been argued on the basis of a fairly reliable effective theory of QCD at high energy, 
that is, a model of color glass condensate (CGC)[1]. 

According to a model of the CGC, the glasma is homogeneous in the longitudinal direction and inhomogeneous 
in the transverse directions. Hence, they form color electric and magnetic flux tubes extending in the longitudinal 
direction. In our previous papers[jj[(|0], we have discussed a possible mechanism for the decay of the color magnetic 
field and have reproduced the instability of the glasma observed in numerical simulations [H, In the discussion, 
only classical gluon's dynamics has been considered. On the othe hand, taking account of quark's dynamics such 
as quark pair creation in color electric field, we have discussed 10] the decay of homogeneous color electric field. In 
particular, we have analytically obtained the temporal evolution of the electric field as well as the number density of 
the quarks. In the analysis, we have not explicitly used the wavefunctions of the quarks in the electric field. Instead, 
we have simply used the formula of chiral anomaly, which already involves the quantum effects of pair creation and 
annihilation. 

In the present paper we generalize our previous discussion [Tpj to electric field which forms an axial symmetric flux 
tube infinitely extending in the longitudinal direction denoted by the coordinate z. We discuss the decay of the 
field which loses its energy owing to the acceleration of produced fermions. These phenomena occur in cylindrically 
symmetric way under the effect of longitudinal strong magnetic field B ^> E. Therefore, we need to discuss an 
azimuthal magnetic field generated by the current of the fermions. Indeed, we can analytically obtain the spatial and 
temporal behaviors of the azimuthal magnetic field and the electric field after the axial symmetric electric flux tube 
is switched on. The fields depend on a cylindrical coordinate r = \/ x 1 + y 2 as well as time coordinate t. As far as 
we know, it is first to obtain the analytic formulation of the azimuthal magnetic field arising in the axial symmetric 
electric field. 

Here we mainly discuss the pair creation of electrons and positrons for simplicity. The generalization to QCD is 
straightforward. Before discussing the pair creation in the electric flux tube, we first give a brief review of the pair 
creation in homogeneous electric and magnetic fields. We assume that the fermions are massless and that the fields B 
and E are spatially homogeneous and parallel (or antiparallel) with each other. They are oriented in the z direction. 
Then, the energies of electrons with charge — e < and positrons with charge e > under the magnetic field B = \B\ 
are given by 



E N = y/p 2 z + 2NeB (parallel) and E N = y/p 2 z + eB(2N + 1) (antiparallel), 



(1) 
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where p z denotes momentum in the z direction and integer N > does Landau level. The term of "parallel" 
( "antiparallel" ) implies magnetic moment parallel (antiparallel) to B. The magnetic moment of electrons (positrons) 
is antiparallel (parallel) to their spin. Thus, electrons (positrons) with spin antiparallel (parallel) to B can have zero 
energy states in the lowest Landau level; their energy spectrum is given by En=q = \Pz\ > 0. On the other hand, the 
other states cannot be zero energy states; their energy spectra are given by \/p 2 + eBM > eBM > eB with positive 
integer M. In a sense they are "massive" states whose masses increase with B. 

When the electric field is switched on, the pair creation of electrons and positrons with the energies En occurs. 
However, it is not probable that any states are produced with an equal production rate. Indeed, high energy states 
are hard to be produced while low energy states are easy to be produced. Thus, the "massive" states are hard to 
be produced by weak electric field E <C B. In particular, they cannot be produced in the limit of B — > oo. Only 
products in the limit are the pairs of electrons and positrons in zero energy states ( £W=o = \p z \ = ). 

Note that owing to the Fermi statistics only fermions with p z = are produced, since the other states with p z 
have already been occupied. After the production of the fermions with p z = 0, their momenta increase with time, 
p z (t) = ±e J dt'E, owing to the acceleration by the electric field. Hence we obtain the momentum distribution 
h(p z ) oc 9(pF(t) —p z )0{p z ) for positrons with charge e > and n(p z ) oc 8(pF(t) +p z )0{—p z ) for electrons with charge 
— e < 0. The Fermi momentum pp is given by pp(t) = e Jq dt'E(t'). Here we have assumed that the electric field is 
parallel to B and is switched on at t = 0. (When E is antiparallel to B, the momentum of the particles increases in 
the different direction from the case of E being parallel B.) 

We note that electrons move to the direction antiparallel to E while positrons move to the direction parallel to E. 
Therefore, both electrons and positrons created by the electric field have right handed helicity when E being parallel 
B, while they have left handed helicity when E being antiparallel B. 

In our discussion of the pair creation, the key point is to use the equation of chiral anomaly, 

e 2 

d t {n R -n L ) = —E(t)B (2) 

where tir ( ) denotes number density of right ( left ) handed chiral fermions; nn t L = (^7o(l ± 75)^)/2 in which 
the expectation value is taken by using a state of electrons and positrons created in the electric field. Here, we have 
used spatial homogeneity of the chiral current j'5, that is, divjs = 0. 

In accordance with the chiral anomaly, the rate of chirality change is given by the product of E and B. It is important 
to note that only particles with right ( left ) handed helicity are produced in the limit B ^> E when E is parallel to 
B (E being antiparallel to B.) Since the number density n of electrons is the same as that of positrons, tlr — 2n and 
nt = when E being parallel to B, while ur = and m, = 2n when E being antiparallel to B. Thus, we find that 
the temporal evolution of the number density is determined by the anomaly equation 20 t n(i) = e 2 \E(t)\B / An 2 . 

After the production of the particles, the electric field gradually loses its energy because it accelerates the particles. 
The energy of the electric field is transformed into the energies e of the particles, so that the energies of the particles 
increase, 

dt (e{t)+ l -E{tf) =0 (3) 

where we have neglected the contribution of magnetic field induced by the electric current of electrons and positrons. 
Later we take into account the contribution when we discuss pair production in an electric flux tube. 

The energy density e(t) of the particles are given such as e(t) = 2n(t)pF(t)/2 = n{t)pp(t). This originates from the 
fact that the momentum distribution h(p z ) ( n(t) = ± J ±o ° dp z h(p z ) ) is given by n(p z ) = n^9(pp{t) — \p z \)Q{\p z \) as 
we have discussed. For example, the energy density of electrons is given such that J_ dp z h(p z )\p z \ = npp/2. 

Now we have three equations to solve for obtaining the pair production rate dtn(t), etc., 

2d t n{t) = ^E(t)B, d t (e(t) + ^E(t) 2 ) = 0, and e(t) = n(t)p F (t) (4) 

with p.p(i) = J Q dt'eE(t'). It is easy to solve the eq(U]) with initial conditions E(t = 0) = E^ > and n(t = 0) = by 
assuming magnetic field B independent on t, 
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FIG. 1: electric field E(t) (solid line) and number density 
n(t) (dashed line) 



E(t) = E cos( 




aeB 



t) and n(t) 



aE Q B\wa{\pf-t)\ 



(5) 




with a = e 2 /47r, where we have taken account of the number density being positive semi-definite. We find that the 



The formula holds rigorously in the limit JS»£o since the particles with masses on the order of B are suppressed 
in the production. In Fig. 1 we have shown the behaviors E(t) = Eocos(7rt/2t c ) and n(t) = 2Eoa ^ tc | sin(7rt/2t c )| 
with Eq = 1, t c = 1 and a ~ 1/137. When the electric field is switched on at t = 0, the pair creation begins to 
occur, so that the number density of electrons and positrons increases. Owing to the acceleration of the particles by 
the electric field, the energy of the electric field gradually decreases and eventually vanishes at t = t c . At the same 
time the number density takes the maximum value. Then, the direction of the electric field is flipped and it becomes 
strong with time. On the other hand, the number density decreases after t = t c . This decrease is caused by the pair 
annihilation of electrons and positrons. Since the direction of the electric field is flipped after t c , the direction of the 
particle acceleration is also flipped. Consequently, electrons and positrons are moved to overlap with each other so 
that the pair annihilation may occur to make the number density decrease. Eventually, n vanishes at t = 2t c . Then, 
the pair creation of left-handed particles begins to occur since E is anti-parallel to B after t — 2t c . Therefore, the 
oscillation of E(t) and n(t) arises. 

We wish to make a comment that although the formula of the chiral anomaly involves all of the quantum effects, 
the mechanism why the number density increases is not explicit in the formula. However, we understand that the 
increase or decrease of the number density originates with the pair creation or annihilation. This is because the results 
shown in this paper have agreed with those obtained by calculations [lj explicitly including the quantum effects of the 
pair creation or annihilation. 

Now we discuss the pair production in an axial symmetric flux tube of electric field, namely the field extending 
infinitely in the longitudinal direction, but has finite width R in the transverse directions. In this case the quantities 
of E, n and pp depend on the cylindrical coordinate r and time coordinate t. ( Since the wavefunctions in the 
Lowest Landau level are given by z m exp(— eB\z\ 2 /A) with z — x + iy and integer m > 0, the charged particles can be 
located within the flux tube for sufficiently large B such as R 1/y/eB. ) Thus, the formula of the chiral anomaly 
takes the same form as the one in eq(j4]), although the electric field E(r,t) and the number density n(r,t) depend 
on r and t. ( Owing to the axial symmetry and homogeneity in the longitudinal direction, we may take divj 5 = 0. 
) Similarly, the energy density e takes the form e(r,t) — n(r,t)pp(r,t) as eqQ, in which the Fermi momentum is 
given by Pf{t, t) = J Q dt'eE(r,t'). We also take into account Maxwell equations such as dtBe(r,t) — d r E(r,t) and 
dtE(r, t) — dr(r.Bfl(r,t)) _ j( r; where J(r, t) represents the current of electrons and positrons flowing in the z direction. 
Bg(r,t) denotes azimuthal magnetic field generated by the current. Therefore, we have the following equations to 
solve, 



life time of the electric field, which may be defined as E(t c ) — 0, is given by t c = f (y ^f— 




2d t n(r, t) = -^E(r, t)B, e(i) - n(t)p F (t), 



(6) 
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and Maxwell equations, 

d t B 9 (r,t) = d r E(r,t), d t E(r,t)= drirBeir ' t)] -J(r,t). (7) 

r 

In order to solve the equations we must represent the current J(r, t) in terms of n(r, t) and pf(t, t). It can be done 
by imposing the condition of energy conservation, 



g f d3 ^ E (r,t) + B g (r,t) + £(r ^ = J d 3 x (E(r,t)d t E(r,t) + B e (r,t)d t B e (r,t) + d t e(r,t)) = 0. (8) 

The condition implies that the energy of the electric field is transmitted into the energies of the azimuthal magnetic 
field Bg and the produced electrons and positrons. The particles are accelerated by the electric field so that the field 
induces the electric current J, which produces the field Bg. 

Using the Maxwell equations in cqj7|), we rewrite the condition such that 



J d 3 x(E(d t E - dr ^ Be ^ ) + d t e^j = J d 3 x(-JE + d t e) = J d 3 xE(en + ^Bp F - J) = 0, (9) 

where we have performed the partial integration in r. The energy conservation must hold for any initial conditions of 
E(r,t = 0) = Eo(r). Therefore, we obtain 



J(r, t) = en{r, t) + ^Pf{t, t)B = 2en{r 7 1), (10) 

since -^pF(r,t)B = J* dt'eE(r,t')B = en(r,t) with n(r,t = 0) = 0. It states that both of electrons with charge 
density —en and positrons with en constitute the electric current J since the velocity of electrons is —1, while that of 
positrons is +1 in the unit of light velocity c = 1. Obviously, the expression J = 2en is very natural. 

Using J(r,t) in eq fTOl) as well as equations in eq((4]), we derive the equation of motion of the electric field, 



a, 2 £(r,() = (a,' 



r 



e 3 B 



(ii) 



where we have used the formula dfJ(r,t) — j^BE(r,t). 

We solve the equation fTTj) with the initial conditions n(r, t 



0) = 0, Bg(r,t = 0) = and E(r,t = 0) 



Eti exp(— r /R ). The initial conditions implies that the electric field forming flux tube with radius R is switched on 
at t = 0, when any particles and magnetic field Bg are absent. Noting that dtE(r, t = 0) = derived from Maxwell 
equations ([7]) as well as the initial conditions, we first obtain the solution E, and then we derive n and Bg, 



E(r,t) 
n(r, t) 
B e {r,t) 



2 r°° 



EqR 



e 2 E R 2 



kdk cos^v^ 2 + m 2 )J {kr) exp(-k 2 R 2 /4), 



16tt 2 

EqR 



sin(tVfc 2 + m 2 ) , 2n2 ,a 
kdk — \ = - J (fcr) exp(-fc 2 ir /4) 



\/k 2 + m? 
sm 



2 poo 



kdk 



\/k 2 + m 2 



Ji(fcr)exp(-ri?74) 



(12) 

(13) 
(14) 



with 



in' 



■B, where Ji(kr) denotes Bessel functions. We have used the fact that the number density n must be 



positive. Taking the limit R — > oo in eq(H2J)^ j eq (IT4"l) . we find that the solutions in ea ([T2^) ~ ea (IT4l are reduced to the 
previous solutions in eq©, which are obtained under the condition that the electric field is homogeneous, in other 
words, the flux tube with infinite width. 

We have shown the temporal and spatial behaviors of electric field E(r, t) in Fig. 2, number density n(r, t) in Fig. 3, 
and azimuthal magnetic field Bg in Fig. 4, respectively, where the parameters m = 1, EqR 2 /2 = 1, e 2 /8ir 2 — 1 and 
R = 1 have been used; the scale r = 10 in the figures corresponds to R. Small and large dots represent the behaviors at 
t = 0.2 and t = 0.6, respectively. We find that the electric field E(r,t) vanishes faster as r is smaller since E(r,t = 0) 
is stronger as r is smaller. The strong electric field produces the particles so much that the number density of the 



■5 



1.5 



ol.O 
W 

0.5 



M • • » 1 1 

5 10 15 20 25 30 



FIG. 2: electric field £(r) at t = 0.2 (small dots) 
and E(r) at t — 0.6 (large dots) 
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FIG. 3: number density n(r) at t — 0.2 (small dots) 
and n(r) at t — 0.6 (large dots) 
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FIG. 4: azimuthal magnetic field Bg(r) at t = 0.2 
(small dots) and Bg(r) at t — 0.6 (large dots) 



particles is larger as r is smaller. We also find from the solutions (fT2|) that the life time of the electric field with finite 
R is shorter than that of the field with R = oo. ( The life time may be defined roughly as E(r — R,t = t c ) = 0, 
although the field oscillates with time. ) This is because typical momentum dominating in the k integration is on the 
order of R . Thus, the life time is approximately given by t c ~ {\Jm? + i? -2 ) -1 . 

As we can see, the quantities such as E(r,t), n(r,t), and Bg(r,t) oscillate with time. In particular, electric field 
E(r,t) never disappear by completely losing its energy. This originates with the fact that the produced particles are 
almost free and never interact with external so that they never lose their energies and momenta. Indeed, we have used 
the momentum distribution of the free fermions. On the other hand, when we take into account their interactions 
with external heat bath and use fermionic momentum distributions at finite temperature, the electric field would 
decay and disappear because the current J would vanish owing to the dissipation of the particle's momenta in the 
heat bath. 

Distinguished point in our analysis is that simply when we adopt an appropriate momentum distribution including 
the effects of particle interactions, we can discuss the disappearance of electric field using the chiral anomaly in the 
limit B ^> E. We do not need to explicitly obtain the wavefuncions of the interacting particles. For this simplicity, 
we are able to explore a mechanism for the disappearance of the electric field. 

It is interesting to see from eq([TT|) that the electric field E parallel to the background magnetic field B possesses an 
effective mass m = aeB/ir with a = e 2 /Air owing to the pair creation of electrons and positrons. This result holds 
only for B E and the pair production of massless fermions. When the fermions have a mass ( indeed, electrons 
have small mass ), we need to modify the equation of the chiral anomaly and momentum distributions of the fermions. 
However, the modification of our results would be minor when the mass is much less than B. We will discuss the 
effect of fcrmion mass on the pair creation in future publications. 

Finally we wish to make a comment on the decay of the glasma produced in high-energy heavy- ion collisions. The 
color electric and magnetic fields generated in the early stage of the collisions form flux tubes. Although we have 
assumed that magnetic field is not a tube-like but homogeneous, we may apply our results to the analysis of the 
glasma decay, in particular, the decay of the color electric field. Then, as we found, the life time t c of electric field 
is approximately given by (V R~ 2 + m 2 ) -1 . In the glasma, m 2 is of order a s gB where gB ~ Q 2 S with saturation 
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momentum Q s ( a s = g 2 /4-7r ~ l/4ir with gauge coupling constant g ~ 1 ), while the width of the flux tube R is 
given by Qj 1 - Hence, we find that the color electric field decays rapidly such as t c ~ Qj 1 ~ 0.1 fm/c ~ 0.2 fm/c for 
Q s = 1 GeV ~ 2 GeV. The smallness of the width is the main origin of the rapid decay. Although this is the result 
obtained in the limit of B ^> E, the result indicates that the color electric field decays sufficiently fast as required 
phenomenologicallydlj. In addition, the result is consistent with our assumption of the independence of B on time. 
Indeed, B decays very slowly [§, [l2T] and does not vary so much at least within the period of Qj 1 . 

To summarize, using the chiral anomaly we have discussed the pair creation of electrons and positrons in an axial 
symmetric electric flux tube under the effect of strong magnetic field. We have analytically obtained the spatial and 
temporal behaviors of the electric field and azimuthal magnetic field induced by the current of the fermions. Our 
calculations have been performed without addressing explicit forms of particle's wavefunctions. However, our results 
hold only in the limit of B 3> E, in which all contributions of "massive" states vanish. Applying our results to the 
glasma produced in high energy heavy-ion collisions, we have found that the color electric field decay sufficiently fast 
to be consistent with QGP phenomenology. 

The author express thanks to Drs. H. Fujii of University of Tokyo and Dr. K. Itakura of KEK for their useful 
discussion and comments. 



[1] N. Tanji, hep-ph/0810.4429; see the references therein. 
[2] J. Schwinger, Phys. Rev. 82 (1951) 664. 

[3] A. Casher, H. Neuberger and S. Nussinov, Phys. Rev. D 20 (1979) 179. 

K. Kajantie and T. Matsui, Phys. Lett. 164 B (1985) 373. 

M. Gyulassy and A. Iwazaki, Phys. Lett. 165 B (1985) 157. 
[4] E. Iancu, A. Leonidov and L. McLerran, |hep- ph/0202270 

E. Iancu and R. Venugopalan, hep-ph/0303204 

For a brief review, see K. Itakura, Prog. Theor. Phys. Suppl. 168 (2007) 295. 
[5] A. Iwazaki, Phys. Rev. C 77 (2008) 034907; hep-ph/0803.0188, to be published in Prog. Theor. Phys. 
[6] H. Fujii and K. Itakura, Nucl. Phys. A 809 (2008) 88. 
[7] H. Fujii, K. Itakura and A. Iwazaki, hep-ph/0903.2930. 

[8] P. Romatschke and R. Venugopalan, Phys. Rev. Lett. 96 (2006) 062302; Phys. Rev. D 74 (2006) 045011. 
[9] J. Berges, S. Scheffler and D. Sexty, Phys. Rev. D 77 (2008) 034504. 
[10] A. Iwazaki, hep-ph/0904.1449. 

[11] T. Hirano and Y. Nara, Nucl. Phys. A 743 (2004) 305; J. Phys. G 30 (2004) S1139. 
[12] T. Lappi and L. McLerran, Nucl. Phys. A 772 (2006) 200. 



